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Abstract 

The effect of spin-charge separation on the kinetic energy of the two-dimensional t-J model is exam- 
ined. Using a sum rule, we derive an exact expression for the lowest possible KE (Ebound) for any state 
without doubly occupied sites. The kinetic energies of the relevant slave-boson and Schwinger-boson 
mean-field states are found to be considerably larger than Ebound- The MF states exhibit complete 
spin-charge separation, and form the basis of a number of microscopic theories. An examination of the 
momentum distribution reveals that the large increse in KE of the MF states is due to excessive deple- 
tion of electrons from the bottom of the band (Schwinger boson states) and holes from the top of the 
band (slave boson states). To see whether the excess KE is simply due to the poor treatment of the 
local constraints, we solve the constraint problem analytically for the Schwinger boson MF states in the 
J — limit. This restores gauge invariance, incorrectly violated in the MF theories. The resulting state 
is a generalization of the Hartree-Fock state in the Hubbard model, but one that includes spin-wave 
excitations, removing a deficiency of the simple HF theory. Even after the constraints are imposed cor- 
rectly, the MF kinetic energy is found to be much larger than Etound- These results support the notion, 
advanced in earlier papers, that spin-charge separation in the MF state costs excessive kinetic energy, 
and makes the state unstable toward recombination processes which lead to superconductivity in d = 2 
and a Fermi liquid state in higher dimensions. 

PACS: 71. fO Fd, 74.25 -q, 74.72 -h, 74.62 Dh 

I. INTRODUCTION 

The possiblity of spin-charge separation in large-C/ Hubbard model and the equivalent t-J model is of 
considerable interest because of its relevance to high-T c superconductivity [1]. Such a behavior is seen in 
d = 1, where the metallic state is characterized by separate spin and charge excitations, not electron- like 
quasiparticles [2]. The metallic state evolves out of the insulating state at half filling, which has no magnetic 
order. By contrast, for d > 2, the ground state is antiferromagnetically ordered at half filling. Therefore, if 
the 2-d Hubbard model were to describe the cuprates, its low-energy behavior would be governed by several 
distinct fixed points: an AF insulator; a metallic, presumably spin-charge separated state without magnetic 
order; possibly a superconducting state; and possibly even a Fermi liquid at lower electron densities. 

A number of microscopic theories are based on the slave-boson (SLB) [3, 4] and the Schwinger-boson 
(SWB)[5] representations of the 2-d t-J model in which the electron field is decomposed into a spin (spinon) 
and a charge (holon) component, subject to a local constraint. Each theory is centered on a mean-field 
state which is supposed to take care of short-range or high-energy processes, leaving behind weak residual 
interactions. The MF state is characterized by independent spinon and holon excitations, and thus exhibits 
complete spin-charge separation. The results have been varied, but not too satisfactory. That none of the 
MF theories succeed suggests that there is a generic reason for the failure - some important bit of physics is 
missing in all of them. Identifying it and understanding the energetics behind the failure is of great interest. 

One obvious source of error, to be considered later, is the local constraint, which is treated on the 
average in all MF theories. We have argued previously on continuity grounds that there is another, and 
more fundamental, reason for the failure [6, 7]. The MF states exhibit maximum spin-charge separation 
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in all dimensions and at all densities. This can not be true, since in high dimensions (e.g., in d = 3), one 
expects a Fermi liquid state to emerge. In d = 1, strong spin-charge separation occurs because electrons 
of opposite spins can not pass each other, and thus are localized as "spins" or "local moments". This 
costs extra kinetic energy. In higher dimensions, electrons can avoid double occupancy by going around each 
other and delocalize, and thereby reduce the kinetic energy. Therefore, spin-charge separation should weaken 
with increasing dimensionality and decreasing electron density. MF states fail because they show too much 
spin-charge separation which costs excessive kinetic energy. They would be unstable toward spin-charge 
recombination into physical electrons, the latter would try to form a Fermi liquid. Recombination will help 
destroy magnetic order because, as the electrons delocalize, the size of the moment is reduced. It will be 
more effective at higher dimensions and far from half filling since more pathways are available, which would 
explain the emergence of the Fermi liquid state in d = 3. 

Spin-charge recombination is a collective process which appears through the kinetic energy term. It has 
been analyzed in d = 2 using RPA in the Schwinger-boson representation. The electron is a "collective" 
excitation. Recombination is strong enough to destroy magnetic order, but the normal state remains weakly 
spin-charge separated [6]. Additional spin-charge recombination, in which a pair of holons combine with a 
pre-paired spin singlets then leads to d- wave superconductivity [7] . In this picture the destruction of magnetic 
LRO, the appearance of superconductivity and the emergence of a Fermi liquid state at high dimensions are 
all caused by spin-charge recombination, as the system evolves away from half filling. 

In this paper we present two exact results which allows us to get more insight into the energetics of spin- 
charge separation. Specifically, we examine the hypothesis that MF states have excessive kinetic energy. 
This issue is quite subtle. How does one know how much KE is too much? With increasing U, kinetic energy 
of the actual ground state will increase, as electrons arc excited out of the Fermi sea. As U — > oo, there is no 
double occupancy, and we expect KE to be quite large (compared with the noninteracting case, Ef ree ^ 7 even 
if the state is a Fermi liquid. By using a sum rule we determine the lowest possible kinetic energy (Eb OU nd) 
for any state of the t-J model, i.e., for any state in which there is no double occupancy. As expected, we 
find that Ef, und >> Ef ree . Interestingly, near half filling, the Schwinger-boson state has a lower KE than 
the slave-boson state. However, as we shall see, the kinetic energy of either state is much larger than Eb oun d- 
An examination of the momentum distribution reveals that excess KE is due to excessive depletion of single 
particle states near the band edges - a situation that is quite likely to give rise to strong fluctuations. 

The second issue involves the largely unsolved problem of the constraints. Treating them on the average 
would be valid if fluctuations were small enough to be treated perturbatively. However, a perturbation theory 
is not likely to converge, and may make matters worse. The constraints imply a local gauge symmetry which, 
according to Elitzur's theorem [8], can not be spontaneously broken. MF theories break this symmetry so 
that elementary excitations, i.e., spinons and holons, are not gauge invariant and therefore not observable. 
The true ground state and elementary excitations are of course gauge invariant. Since gauge invariance can 
not be restored perturbatively it is hard to judge the correctness of the theory, even at a qualitative level. 

One way to avoid the constraint problem is to consider the Hartree-Fock (HF) states of the Hubbard 
model [9], in which constraints are treated exactly. These states are related to the Schwinger boson MF 
states. But the HF theory has its own shortcomings. The spin-wave excitations are missing at the HF level, 
so that low-energy behavior is not described correctly. For example, magnetic order persists all the way upto 
a T c ~ U in all dimension, violating the Marmin- Wagner theorem. On the other hand, for the Heisenberg 
model, the Schwinger boson MF theory is qualitatively similar to the spin-wave theory, which is based on 
the Holstein-Primakoff (HP) representation and is gauge invariant [10, 11]. But HP theory applies only to 
half-filling, there is no charge excitations. 

Here we solve these problems in the U = oo limit. This limit is of interest because the mapping onto 
the t-J model is exact, so a direct comparison between SWB and HF results can be made. By integrating 
out the constraints in the Schwinger boson theory, we obtain a generalized HP representation which is valid 
away from half filling. This, in turn, leads to a generalized Hartree-Fock Hamiltonian, which include both 
charge and spin excitations. The procedure restores gauge invariance, and the inclusion of the spin-waves 
takes care of the deficiencies of the HF theory. We find that even after the constraints are taken into account, 
the nature of the state is not changed, nor is the degree of spin-charge separation. Kinetic energy is reduced, 
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but continues to be excessive compared with Ebound, which supports our hypothesis. 



II. THE MODEL AND SPIN-CHARGE SEPARATION 



We start with a discussion of the meaning of spin-charge separation. The Hubbard model is characterized 
by nearest-neighbor hopping t, and on-site repulsion U. For large U, doubly occupied sites can be projected 
out, which leads to the t-J Hamiltonian: 



where the summation is over nearest neighbors. The second term describes superexchange, with J = 4t 2 /U. 
The first term describes hopping, and has the same form as that in the Hubbard model. However, c| CT creates 
a "projected electron", i.e., an electron of spin a, only if the site i is unoccupied. At half-filling, one needs 
to consider only the exchange term. 

In d = 1, the Hubbard model has been diagonalized by Lieb and Wu [12], using the Bcthe Ansatz. The 
occurrence of spin-charge separation can be seen most clearly in the large U limit where, as shown by Ogata 
and Shiba [13], the ground-state wave function for N e electrons can be approximately written as: 



Here ip c harge(xi,X2, ■•) is a dctcrminantal wave function for N e spinless fermions, where Xi are the electron 
coordinates, irrespective of spin. And -0spm is the Bethe's solution for the Heisenberg model for N e spins (with 
the holes removed). The wave function at smaller U is more complicated, but is continuously connected to 
the large U solution in the renormalization group sense. Of course, the low-lying excited states also similarly 
decompose - i.e., the system is characterized by independent spin and charge excitations. A good account 
of the meaning of Eq. (2) and its consequences is given by Anderson [2] . 

Usually spin-charge separation is associated with the Luttinger liquid state, which is a metal and has no 
magnetic order. However, Eq. (2) clearly suggests that the phenomenon is more general. Any wave function 
which has the form of (2), or is continuously connected to such a wave functions, exhibits spin-charge 
separation. Consider the following cases. 

a) Half-filled case: If we continue the Lieb-Wu wave functions to half-filling, the form (2) will clearly be 
preserved. The (spinless) charge band becomes filled, and the system turns into an insulator. A moment's 
reflection shows that such spin-charge separation occurs at half filling in all dimensions. The reduction of the 
large-U Hubbard model to the Heisenberg model is precisely due to spin-charge decoupling. There is exactly 
one charge per site, which can be viewed as a filled band of spinless fermions. For d > 2, the insulating state 
is magnetically ordered. 

b) Mean-Field States: Similarly, mean-field wave functions for the t-J model also exhibits spin-charge 
separation. In this case the electron operator is represented as: 



The spin-operators are then represented as: = bj^bii, and Sf — \{b\^b^ — b\^>i{). In the Schwinger 
boson representation, the spinon is a boson and holon, a fermion. The opposite is true in the slave boson 
representation. In cither case, constraints are treated on the average, and a MF approximation leads to the 
Hamiltonian: 




(1) 



1p = ^charge^. 



spin • 



(2) 



4 = 4a. 

where b\ a creates a spinon and hi creates a holon, subject to the constraint 

&! T &i T +&t i & ii +/i}/i i = l. 



Hmf = Hi 



spinon 



+ Hhol, 



(3) 
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Note that all eigenstates of Hmf are of the form of Eq.(2). This is also true in the presence of arbitrary 
spinon-spinon and holon-holon interactions. 

In general bosons condense, leading to long-range order (LRO). In the Schwinger-boson case, the ordering 
is magnetic. The mean-field theory gives quite a good description of the Neel state at half filling [10]. Away 
from half filling, the AF insulator evolves continuously into a metallic state with a spiral LRO, and eventually 
to a ferromagnetic metal at small J/t [5]. In the metallic region, ipcharge in the SWB theory is again a spinless 
fermion determinant. 

c) Hartree-Fock states: The magnetically ordered states obtained by the Hartree-Fock treatment of the 
Hubbard model [9] constitute a third example of spin-charge separation. The HF decomposion of the U term 
leads to the quadratic Hamiltonian 

Hhf = H -2UY i im.S i , (4) 

i 

where Ho is the hopping part and m, =< Sj > is the local magnetization, which is determined self con- 
sistently. To describe spiral states, it is convenient to break the spin rotational symmetry in the xy plane 
and choose < >= me"' T as the order parameter, where Q is the spiral wavevector which continuously 
connects the AF state (Q = (n, it)) at half-filling, to the ferromagnetic state (Q = 0) which occurs at large 
US/t. The theory has been described in detail elsewhere [9]. The HF approximation leads to two magnetic 
bands separated by a gap of size 2Um. At large U, low-energy states are confined to the lower band. How- 
ever, there is no spin-degeneracy; each one-particle state can accomodate one electron, not two. Therefore, 
low-energy wave functions are one-component (i.e., spinless) determinantal wave functions, similar to ipeharge 
in Eq. (2). The HF states are spin-charge separated because they are closely related to the corresponding 
states of the Schwinger boson MF theory. 

To see this, consider H HF for Um >> t. Then, the spin of an electron will locally be parallel to the local 
moment mi. Let // and g\ create electrons with spin parallel and antiparallel to m ; , respectively. The bare 
electron operators, which are quantized along the z direction, can be expanded as: 

<Ha = u itJ f ia + sgn{a)u*_ a g itJ , (5) 

where cr = 1 for f and —1 for j. The complex numbers u ia depend on the direction of mj, and satisfy the 
normalization condition 

5^10^ = 1, (6) 

cr 

Suppose, irij = (m, 0, </>), then a simple representation is: 

u t] = cos^/2), Uil = sin(fV2)e^\ 
Note that the second term in Hhf is diagonal in the /, g representation: 

-Um^fUi-gUi)- 

i 

For U — > oo, we can neglect the g terms which constitute the upper band. The similarity with the Schwinger 
boson representation is obvious, with u ia playing the role of bi a and fa the charge part of the electron 
operator. The holon operator is obtained by a particle- hole transformation: fa — > h\. Eq. (6) now plays 
the role of the constraints. The difference is that Vs are operators, where as u's are c-numbers. However, 
in the MF approximation, 6's condense, and we can take u ia «< b ia > . For U = oo, the ground state 
is a ferromagnetic metal away from half filling. Then = m, (and hence ui a — u a ), is independent 
of i. The HF wave function, written in terms of the original electron coordinates, (those created by bare 
operators a ia ), spin-charge separates, with ip c harge again given by the one component Slater determinant, 
and is independent of m. ipspin describes the condensed part and is just a constant. 

For finite but large U, there will be small admixture of <?'s via the hopping term which can be eliminated 
perturbatively, to obtain contribution of order J — the exchange contribution. However, as in the t-J model, 
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such corrections are due to spin-spin interaction and thus does not change our conclusion about spin-charge 
separation. There are of course some important differences between the two theories. In the SWB MF 
theory, constraints are treated on the average, whereas in HF theory they are treated exactly. Similarly, in 
the HF theory, the low energy spin excitations are absent, and hence physics of the singlets is not included. 

Spin- Charge Separation and Local Moments 

At half- filling, it is natural to associate spin-charge separation with the formation of local moments. In 
d = 1 , evidently the local moment picture persists away from half-filling since ipspin continues to be the wave 
function of a "squeezed" Heisenberg antiferromagnet [13, 2]. At large U, two electrons of opposite spin can 
not go past each other. Hence, electrons are localized within a region of size (1 — <5) _1 , and the size of the 
moment is approximately (1 — S)/2. On the other hand, if they can delocalize there will be no moment, 
nor will there be spin-charge separation. Therefore we can associate spin-charge separation with "localized" 
moments. Evidently this picture is consistent with Schwinger boson MF states, as well as HF states. Note 
that spin-charge separation is associated with the amplitude, and docs not require the existence of LRO 
(phase coherence) . The latter disappears at any finite T in d = 2 (Mermin- Wagner theorem) . But the wave- 
function remains spin-charge separated since, at the MF level, the moment vanishes only at T ~ U. Because 
moment formation involves a type of localization, it would cost kinetic energy. Interestingly, from this point 
of view a transition from a Luttingcr liquid to a Fermi liquid can be considered a form of derealization 
transition. 

III. SPECTRAL SUM-RULE AND LOWER BOUND ON KINETIC ENERGY 

In this section we obtain an exact expression for the lowest possible KE in the t-J model. It is useful to 
express the hopping part of the t-J Hamiltonian in the momentum space 

H hop = ^e(k)4 CT c kCT , (7) 

k,<T 

where e(k) = — 2t(cosk x + cosk y ) is the free-electron hopping energy, and c\. is the projected electron 
operator in the momentum space. The average hopping energy (per site) is then given by 

£fcm = ^^e(k)n c (k), (8) 

k 

where n c (k) = n ccr (k) is the (total) electron momentum distribution, with n ca (k) =< cj^Cka >• Eq.(8) 
is exact. Hence any errors in Eun are due to approximations used in calculating n c (k). The kinetic energy 
of the Hubbard model is given by the same formula, except that the corresponding momentum distribution 
differs from n c (k) by terms of order J, and J 5. Here, we focus on the kinetic energy as defined above. Of 
course, for U = oo (i.e., J = 0), the mapping is exact. 

Usually one is interested in the singularities of n c (k) near the Fermi surface. However, in the present 
case, we need to know n c (k) for all k. In order to minimize the kinetic energy, the system will try to make 
n c (k) as large as possible near the bottom of the noninteracting band e(k), and as small as possible near 
the top. Hence states with k near the band edges are quite important. This has serious implications for 
the mean-field states. The assumption underlying the microscopic theories is that high energy processes, 
i.e., those involving electrons with k far from the Fermi sea, are accurately taken care of at the mean-field 
level, leaving behind residual interactions which describes low-energy physics. Therefore, if a MF state does 
not yield accurate values of n c (k) near the band edges, it is likely to be unstable toward fluctuations which 
transfer electrons from the top to the bottom of the band, and thereby lower the kinetic energy considerably. 

A. Spectral Sum Rule: 

Now, the noninteracting state has the the lowest KE, with n c (k) = 2 below the Fermi level ep. However, 
this is not a good reference state since any interaction would increase the kinetic energy by exciting electrons 
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out of the Fermi sea and, as we see below, at large U, a substantial part of of the increase comes from 
avoiding double occupancy, which may have nothing to do with spin-charge separation. Fortunately, we 
can determine the lowest possible kinetic energy in the t-J model. Since the latter does not have double 
occupancy, the corresponding increase in KE is accounted for at the outset. 
Formally, we have 

f dz 

rw(k)= / — A CCT (k,z)/(z), (9) 

where A CIJ {k,uj) is the spectral function for projected electrons which obeys a sum rule [14]. Here we give 
a derivation which is useful for this paper. Let us define the momentum distribution for physical holes (not 
to be confused with holons): p ccr (k) =< c^c^ >. This is obtained from 

P CCT (k) = J ^A ca (k,z)(l - f(z)). (10) 
Hence, total spectral weight for spin a is given by 

For usual fermions operators, the anticommutator equals unity, giving the usual spectral sum rule. This is 
not true for projected electrons. To evaluate the right hand side, we use the representation of the electron 
operator in terms of spinon and holon operators which, in momentum space, gives 

q 

Then the average value of the anticommutator is given by 

< {c kff , clj >= — J2(< b\ a b ia > + < h\hi >) = n ba + 5, (12) 

i 

where is the average number spinons (per site) of spin <r, and 5 is the average number holons per site. 
It is useful to separate into a charge and a spin part. Let = rib/ 2 + °rn, where m = (n^ — ^6j.)/2 is the 
average magnetization, and rib = n ba = 1 — S (from the constraint condition). Then we have 

n ca (k) +p C(J (k) = 1(1 + 6) +am. (13) 

We stress that these sum rules are exact. Since the right hand side depends on the spinon and holon 
densities the sum rules are valid in cither representation, and in all dimensions, irrespective of whether the 
local constraint is imposed exactly or on the average. Summing over a we obtain the total electron spectral 
weight: 

n c (k)+p c (k) = (l + 5). (14) 

The remaining (1 — S) resides in the upper Hubbard band. Note that, in contrast to bare electrons, the 
spectral weight (Eq. 13) depends on charge and spin densities. 

The sum rule, Eq. (14), implies that each momentum state k can accomodate at most 1 + 5 electrons 
(compared with 2 for a noninteracting system). It follows that, at half filling, in order to accomodate all the 
electrons, one must have n c (k) = 1 for all k, which is the correct result. In the Schwinger boson MF theory 
this can be viewed as a filled holon band. The kinetic energy per electron is maximum (= 0), as compared 
with ~ —t, for free electrons. 

B. Lower Bound for the Kinetic Energy: 

Away from half filling, the exact n c (k) is not known. But it is clear that states with low kinetic energy 
will have n c (k) large near the bottom of the band and small near the top. Indeed, the lowest possible kinetic 
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energy for the model is obtained by filling up a modified Fermi sea with exactly n c (k) = (1 + 5), upto a new 
Fermi level cfi. Above cfi, n c (k) = 0. This gives an exact lower bound for the kinetic energy per site. 

E bound = ^ ]T < k )°( e Fi - e ( k ))> ( 15 ) 

k 

where e^i is determined from 

k 

where n c — 1 — 6. The modified Fermi sea is obtained from the noninteracting one by exciting 1 — 6 electrons 
from each state below the (noninteracting) Fermi level ep to a state above. Therefore cfi > ep, and the 
kinetic energy is correspondingly larger (see Fig. 1). This cost in kinetic energy is due to the removal 
double occupancy. In other words, Ebound is the lowest possible kinetic energy for any state in which double 
occupancy is forbidden. 

Apart from the prefactors, the expression for Ebound at a density n c = 1 — 5 is identical to that of a 
noninteratcing system at a higher density n' — 2n c /(l + 6): 

Ebound(n c ) = — Ef ree ( ^ -j-^ ) . (16) 

The volume of the modified Fermi sea is 2/(1 + 5) times the free electron Fermi sea. 

What sort of a state corresponds to E bound ? Since it evidently resembles a free-fermion state with a 
sharp Fermi surface, let us consider the following wave function: 

I* >= ZU k<kF .o. cljvac >, 

where Z is a normalization constant. Since c^ a represents a projected electron, it is easy to see this is the 
Gutzwiller wavefunction in the infinite-?/ limit. Now, the projected electron operators for arbitrary spins 
satisfy the following commutation rules 

{^, c L, (T ,}-^E e4(k '~ k) ' r *[ 6 t' & - + (5 -' /i l^' 

i 

Notice that right-hand side, which arises from the on-site anticommmutators, consists of Fourier components 
of charge and spin-density operators. Suppose, as an approximation, we replace the right-hand side by its 
average value and assume that there is no magnetic or CDW order. Then, we obtain 

{Ckcr,4.',cr'} ~ £<W'<5kk', 

where £ = (1 + 5) /2. In other words, Cko- is effectively replaced by £^ times an usual fermion operator which 
carries spin. If kp is chosen to yield the correct electron density in the same approximation, then the new 
Fermi surface agrees with the modified one, and the energy is given by Ebound- Since spin- and charge-density 
fluctuations are completely suppressed, this state does not have any localized moment, and hence, there is 
no spin-charge separation. 

IV. KINETIC ENERGY OF THE MEAN-FIELD STATES 

The kinetic energy of an actual state is greater than Ebound- Now, the momentum distribution for the 
MF states is given by a convolution: 

n ^ = ^E n "( k+c i)( 1±n ^ c i))' 

q 
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where rib(k) and n/ l (k) are the momentum distribution functions for spinons and holons, respectively. The 
minus (plus) sign corresponds to the Schwinger (slave) boson representation. Consider the MF states at 
T = 0. 

a) Slave-boson State: We consider an RVB type state in which bosonic holons condense at k = 0, giving 

n c (k) = l-<5(l-n 6 (k)). (18) 

Spinons form their own Fermi sea, identical to that for noninteracting electrons, and n& (k) = 2 for occupied 
and for unoccupied states. Hence, electrons have the same Fermi surface, but n c (k) equals 1 + 8 below, 
and 1 — 8 for all k above the Fermi level all the way to the top of the band. Since J2k e l ) = 0; the kinetic 
energy per site is given by 

E slave =8E free (l-8). (19) 

b) Schwinger-boson State: We consider the ferromagnetic state (small J) since the kinetic energy is 
lowest in this case. (At larger J, the kinetic energy, as defined here, is larger due to the competition with 
antiferromagnetic fluctuations). Spinons condense at k = 0, so that 

n c (k) = (l-<f)(l-n/,(-k)). (20) 

The holons form a Fermi surface, centered at the zone corner, with n/j(k) = 1 in the occupied region, and 
zero in the unoccupied region. The electrons have the same Fermi surface, so the electron Fermi energy is 
much higher than the noninteracting one. On the average, there are only 1 — 8 electrons below this Fermi 
level (even less than the number at half filling!) and none above. In this case the kinetic energy per site is 

E sch = ^±E free (2S). (21) 

Fig (1) and (2) shows n c (k) and the kinetic energy for various states. First, note that all states have 
the correct kinetic energy (= zero) at half filling. Away from half filling, KE for either MF state is much 
larger than E^ oun( i 1 by about a factor of two in the small- 8 region. Surprisingly, for 8 < 1/3, the Schwinger 
boson state has a lower kinetic energy than the slave boson state. One expects Schwinger boson theory to do 
better as far as exchange energy is concerned since it works well for the Heisenberg model. It appears that 
it is also better for the KE at low doping, the region of interest for the cuprates. In the slave boson case, 
the Fermi surface coincides with the noninteracting one, and below the Fermi level, n c (k) is the maximum 
(= 1 + 8) allowed by the sum rule. However, above the Fermi level, n c (k) is large (= 1 — 8) all the way to the 
top of the band, whereas, it is zero for the other states. This costs considerable amount of kinetic energy. 
In the Schwinger boson case, the top of the band is empty, which is fine, but the kinetic energy is still large 
compared with Eb OU nd because the Fermi energy is much higher, and because n c (k) = 1 — 8 in the occupied 
region, smaller than value allowed by the sum rule, all the way to the bottom of the band. 

c) Hartree-Fock States in the Hubbard Model: How are these results affected by the constraints? That 
there is an error due to the average treatment of the constraints can be seen from the average electron density 

n c = < >= jf Yl nc ( k )' 

a k 

which does not equal its actual value 1 — 8 in the MF theory. Instead, it equals (1 — S) 2 in the Schwinger boson 
state, and (1 — S 2 ) in the slave boson state. To assess the role of the constraints we consider the HF states of 
the Hubbard model in which constraints are treated exactly. We consider U = oo limit since mapping onto 
the t — J model is exact (with J = 0). The HF ground state is fcromagnetic, with the up-spin band filled up 
to the holon Fermi level. The down-spin band is empty. This is the same ferromagnetic state found in SWB 
theory, and has the same fermi surface. However, now n c (k) = fi-f(k) + «j.(k) — 1 in the occupied region, 
instead of 1 — 8. Since, as shown rigorously in section V., the states are essentially the same, the difference 
between the two results is a measure of the error (due to the constraints) in the Schwinger-boson theory. 
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Since the Fermi levels for the HF and Schwinger boson states coincide, the energies are proportional: 

E ac h = (1 - S)E hub . 

Now, for all the states, Ek in < at finite 6. Hence, Ehub is less than E sc h (and also less than E s i ave ) for 
all 6 (see Fig. 2). At small 6, Ehub scales with 6, therefore the corrections due to constraints are small, i.e., 
of order 6 2 . To put it another way, since E sc h/ Ehub = 1 — 6, the effect of the constraints in the physically 
important small 6 region is not large. 

We see that kinetic energies of the MF and HF states are much larger than Ebound, the lowest possible 
KE without double occupancy. In the small 6 region, E), oun d/ Emf is between 1.5 to 2, i.e., the difference 
is comparable to the energy itself (Fig 3). This does not mean that E^ound corresponds to the actual state 
of the system. The key point is that n c (k) considerably smaller (by an amount oc <5 - the scale of KE) than 
the allowed value near the bottom of the band for the SWB and HF states, and considerably larger near the 
top of the band for the slave boson states. In the actual ground state, n c (k) is sure to depend on e(k), and 
is likely to be close to 1 + 6 near the bottom, and close to zero near the top of the band. This means that 
fluctuations (in the MF states) that transfer electrons from near the top of the band to near the bottom will 
be dominant. 

V. CONSTRAINTS: FROM SCHWINGER BOSONS TO SPIN WAVES 

In this section we solve the constraint problem in the U = oo limit. It is convenient to use the functional 
integral technique. The Bose fields can be expressed as: bi a = r <T e % ^ i,J . The constraint condition is given by 

r ?T + r il + K h i = !> 

which does not depend on the phases. The action is invariant under the gauge transformation: bi„ — > 
bi<je lXi , hi — > hie lXi . We can choose \i = — <AiT' ( or ' equivalently select a gauge: <j>^ = 0). The action is 
now independent of <f>^, which can be integrated out. The constraints are not effected. Next, we introduce 
the constraints into the functional integrals as 6 functions, and use the latter to integrate out r^. This is 
equivalent to the replacement = = (1 — r|i — h*hi) 1/>2 in the action. The remaining degrees of freedom 
define the generalized HP representation. It is easy to see that the corresponding Hamiltonian is obtained 
from the Schwinger boson Hamiltonian by using the following substitutions: 

bii — ► hi, hi — > h i; 

^^^(l-btbi-hjh^i 

The spin operators are represented as 

S+ = (1 - bjb, - bjh^ bj . 
Sf ^(l-hjh^-btb,. 

In the absence of the holes, one recovers the usual HP representation for the Heisenberg ferromagnet. The 
generalization to the antiferromagnetic case is straightforward. As in the usual HP representation, the 
square-root represents an anholonomic constraint. This is of course not as stringent as the constraints in the 
Schwinger-boson (or slave-boson) representation. The electron operators are given by 

4i = blh, 

C f T = (l-bthi-hth^h;. 

In the last expression, we can drop the hole density operator within the square-root since it will not contribute 
(because h is fermionic). This gives, cj T = (1— bjbj)5hj. Since the constraints have been integrated out, there 
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is no gauge freedom left. Note that the b and h are gauge invariant since they are proportional to gauge- 
invariant operators S + and <y , respectively. As in the Heisenberg case, b^ creates a spin- wave excitation - 
a magnon. Similarly, the operator creates a spinless hole excitation (not the unobservable holon of the 
SWB representation). Other difficulties associated with the constraints also disappear. For example, the 
local density of electrons is given by 

which correctly yields, n c = 1 — S, in contrast to, n c = (1 — S) 2 in the Schwinger boson MF theory. 
For J = (U = oo) the t-J Hamiltonian, in the generalized HP representation becomes 

H = tg[{(l - bjbi)^! - btb,-)' + btb^hjh, + h.c.}. (22) 

ij 

Note that the total no of holes (Nh) and spin-wave excitations (Nb) are separately conserved. This follows 
from the conservation of the total electron number, N c = N — Nh, and total magnetization M z — J2i &f = 
r;(N — Nh) — Nb in the parent Hubbard model. We will keep N c fixed, which fixes Nh- But Nb is not fixed 
in the ordered state, but varies from to N c , where Nb = N c corresponds to zero magnetization. 

We can diagonalize the Hamiltonian separately in each Nb subspace. The Nb — subspace only has holes, 
and is of particular interest. In this case, hj — ctp All electrons have spin up, the average magnetization, 
m = (1 — 6)/2, has the largest possible value. Within this subspace, the Hamiltonian reduces to 

H h = t^204^j + h.c), (23) 

ij 

with no residual interaction. Note the this is nothing but the U — oo Hubbard Hamiltonian in the Hartree- 
Fock approximation in the fully spin-polarized subspace. In particular, it contains the ferromagnetic HF 
ground state (the generalized Nagaoka state) considered previously. We have thus shown that HF theory 
can be obtained from the Schwinger boson theory, and it takes the constraints into account exactly. 

All wave functions corresponding to Hh are one-component Slater determinants, describing the charge 
sector. These are also exact wave functions of the full Hamiltonian (actually for all U) since in the Nb = 
subspace, Hh is the exact Hamiltonian. It is not known rigorously whether the Nagaoka state is the true 
ground state at finite S. A possible answer is that the true ground state belongs to a different Nb subspace, 
with magnitization m = |(1 — S) — rib, which decreases continuously to zero 5 cr . Here, n& = Nb/N is the 
magnon density. 

For finite rib, we need to include the magnons. For small rib, we expand the square root to leading order 
in rib, which gives 

H = t ^(hjhj + h.c) +tj^[(b]b i iblb, - ibjb,)hjh, + h.c.}. (24) 

ij ij 

The second term includes spin-wave effects, and extends the treatment of the U = oo Hubbard model beyond 
the Hartree-Fock level. 

Note that there is no free magnon terms in the Hamiltonian. This is very different from the Heisenberg 
ferromagnet, where the lowest order (1/S) term is quadratic in bosonic operators. Here, magnons (spin) 
interact with holes (charge) at the lowest level. The simplest approximation is to do a MF decomposition 
of the quartic term, leading to quadratic Hamiltonians for holes and magnons: H = H b + H h . The hole 
Hamiltonian has the HF form (Eq. 23), but with a renormalized hopping parameter 

t h =t(l + B-n b ), 

where, B =< bjbj >. The bosonic mean-field Hamiltonian is given by 

H b = tD ^[btbj + btbi - bjbi - b]bj], (25) 

ij 
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where D =< hjhj >. The magnon Hamiltonian has the same form as the one in the Heisenberg model. The 
magnon energy is then given by, 

w m (k) = 2tD(cosk x + cosky — 2). (26) 

It is easily shown that D < is negative, so that w m (k) > 0. Therefore, there are no magnons at T = 0, 
and the Nagaoka state found in the HF approximation continues to be the ground state at this level of 
approximation. At finite T, magnons will proliferate just as in Heisenberg ferromagnet, and destroy long- 
range order. This is the correct result (Marmin- Wagner theorem). In principle, Takahashi's modified spin- 
wave method [15] can be used to describe the system at finite T. 

Electron Green's Function and Missing Spectral Weight: 

To obtain the Green's function for the up electrons, we expand the square-root. Now, at the MF level, 

T T < c iT (r)c] T (rO >= T T < (1 - \b\(r)b t (r) ib](r')b,(r')) >< hj^h^r') > . 

The average of the density operators simply gives the magnon density rib, which is uniform. Hence, the 
Green's function is proportional to the hole Green's fumction: 

G cT (k,o;)= . l ~ H * (27) 

ILU + 6h(k) 

where e/ l (k) = 2th (cos k x + cos k y ) — fih is the hole spectrum, which includes the chemical potential fih- This 
term is analogous to the condensate part in the Schwinger boson MF theory [16]. 

The Green's function for down electrons is a convolution of magnon and hole Green's functions, and is 
given by [17] 

U ' iu;-w ro (k + q) + e h (q) 1 ; 

This term is the analog of the non-condensate part in the SWB theory. The corresponding spectral function 
is nonzero only for u> > and is incoherent. It is absent in the HF theory. In the HF case there are two 
bands, described by simple poles. The up-spin band is the lower Hubbard band. The down-spin band is the 
upper Hubbard band which represents the doubly occupied sector amd is separated by U — > oo; it is is absent 
in our case. In the simple HF theory, the spectral weights for each band is unity. This does not agree with 
the sum rule, which gives 1 + 6 and 1 — 6, respectively. The reason for the discrepancy is that magnons are 
missing in the simple HF theory. The expression for Gi (Eq. 28), is precisely this (low- frequency) magnon 
contribution, with a total spectral weight of rife + 6. Together with the spectral weight of G-f, this gives a 
total spectral weight of 1 + 6, in agreement with the sum rule. 

The Green's function in the spin- wave theory is very similar to that in the Schwinger boson theory [16]. 
However, in the latter case, magnetic ordering is along the x direction, so that the Green's function for up 
and down electrons are identical. For a proper comparison we need to take the average: G c = (G c -\ + G c {)/2. 
Although the structure of G c in the two cases is the same, the spectral weight distribution between the 
coherent (pole) and incoherent pieces is different. For a comparison, it is convenient to express these in 
terms of the magnetization m which is the same in both theories, and equals (1 — 6) at T = 0. Recall that 
n c (k), and hence, the kinetic energy, is determined by the negative frequency part, i.e., by the pole piece 
(Gci), which carries a spectral weight = |(1 + 6) + m, which is greater than the (total) weight in the SB 
theory = 2m. In particular, it reduces to the HF value (= l)at T = 0. Therefore, HF theory gives the best 
kinetic energy. 

VI. CONCLUSIONS 

In this paper we have shown that mean-field states of the t-J model exhibit complete spin-charge sepa- 
ration in the sense of Eq. (2), and in d = 2, have high kinetic energies. This supports qualitative arguments, 
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given in previous papers, that spin-charge separation localizes electrons which costs kinetic energy. The 
determination of Ebound is a key result since it includes excess KE which arises simply from the removal of 
double occupancy, but not necessarily from spin-charge separation. The calculation of n c (k) is also inter- 
esting, because it shows that mean-field values differ substantially from what is allowed by the sum rule at 
the band edges. Therefore fluctuations that transfer electrons from the occupied to the unoccupied regions, 
can, in principle, lower the energy significantly, i.e, by an amount comparable to KE itself. 

The sum rule and the expression for Ebound are valid for arbitrary J. The kinetic energies presented here 
are for U = oo ( J = 0) for which KE for MF and HF theories can be compared directly. Extension to finite 
J is not difficult since the MF problem has been solved exactly, but numerically. At finite J, the KE of 
the t-J model would be larger since the two terms compete. But conclusions of this paper are not going to 
change. Of course, J will be important to describe the spin singlets and for superconductivity. 

The second important result is the mapping of the t-J model onto a generalized Holstein-Primakoff model. 
This allowed us to solve the difficult problem of the constraints for the SWB MF state at U = oo. Our 
results show that constraints do not change the nature of the state - it remains fully spin-charge separated 
and has the same magnetic ordering. Extending the method to finite J is somewhat more complcated, but 
can be done. It will be useful to study fluctuations such as spin-charge recombination within the generalized 
HP model so as to take into account the constraints at the outset. 

The mapping also establishes the cannection between the SWB theory and the HF approximation. The 
latter has long been used to study magnetic ordering in strongly correlated systems. However, HF approx- 
imation fails at finite T (or, at finite frequency) because spin-waves are not included. For example, T c at 
which magnetic order vanishes scales as U. The inclusion of the spin-waves via the mapping solves this 
problem, as spin- waves destroy LRO in 2-d at any finite T, and in d = 3 at a temperature ~ t (as U — ► oo). 
It also demonstrates complete spin-charge separation in the HF plus spin-wave theory. This method will be 
useful in studying itinerant ferromagnetism in d = 3. 

Our results suggest that fluctuations that result in spin-charge recombination are likely to be strong. 
Such processes lower KE by moving electrons from above the Fermi level to the bottom of the band. In 
previous papers we have treated recombination by RPA in which the physical electron appears as a fermionic 
collective mode [6]. In d = 2, it helps destroy the magnetic LRO, but, for 5 not too large, the Fermi liquid 
state is not recovered. For finite J, additional spin-charge recombination between a pair of holons with 
preformed pair of spinons leads to d-wave superconductivity [7] . 

For the U = oo case, the dcstabilization of the HF state has been demonstrated by Shastry, Krishna- 
murthy and Anderson [18]. They considered the destruction of ferromagnetic LRO away from half-filling. 
Consider the ferromagnetic ground state with all spin up. Since = 1 in the occupied region, the only 
way kinetic energy can be lowered is by destroying an up spin electron at the Fermi level, and creating a 
down-spin electron at the bottom of the band. But creating a bare down-spin electron would cost U. Instead 
one must therefore create a projected down-spin electron. They showed that the MF state becomes unstable 
with respect to such spin modes. Note that in HP representation, (or SWB representation), the projected 
electron is a collective object created by cj^ = N^ 1 / 2 J^q t>k+ q o-^q- This process therefore corresponds to 
spin-charge binding, and is the same one considered by us. 
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FIGURE CAPTIONS 

Fig.l: 

Momentum distribution for various states as a function of energy. The slave boson case has a high kinetic 
energy because n c is large all the way to the top of the band. For the other three states the top of the band 
is empty, but for the Schwinger boson case the energy is again large because n c is small at the bottom of 
the band. 

Fig.2: 

Kinetic energy as a function of hole density d. The HF energy is lower than either Schwinger boson 
or slave-boson case. However, Eb oun d is much lower than Ebound- The discrepancy between the HF and 
Schwinger boson energies measures the error due to poor treatment of constraints. 

Fig.3: 

Ebound/ Ehub as a function of hole density. At small 5, the HF energy is larger by about a factor of two. 
(Note both energies are negative). 
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